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We derive an analytical expression for the longitudinal magnetoconductivity cr zz in layered con- 
ductors in presence of a quantizing magnetic field perpendicular to the layers and for short-range 
in-plane impurity scattering in frame of the quantum transport theory. Our derivation points out 
quite unusual temperature and magnetic field dependences for Shubnikov-de Haas oscillations in the 
two-dimensional limit, i.e. huo c 3> *int, where t is the interlayer hopping integral for electrons, and ui c 
the cyclotron frequency. In particular, when hui c 3> 4irt and huj c > 27rF M (here T M is the value of the 
imaginary part of the impurity self-energy at the chemical potential p), a pseudo-gap centered on 
integer values of p/hui c appears in the zero-temperature magnetoconductivity function a zz (p/hui c ). 
At low temperatures, this high-field regime is characterized by a thermally activated behavior of 
the conductivity minima (when chemical potential p lies between Landau levels) in correspondence 
Q [ with the recent observation in the organic conductor /3"-(BEDT-TTF) 2 SF 5 CH 2 CF2S03. 
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t In the past few years, numerous studies of the effects of a quantizing magnetic field in low-dimensional systems, 
C^) i especially in quasi two-dimensional (2D) organic conductors have been performed (see e.g. [1] and references therein). 
The quasi-2D conductors consist typically of highly conducting planes with a small interplane quasiparticle dispersion. 
In presence of a magnetic field H perpendicular to the layers, the motion in the conducting planes is quantized. The 
quasiparticule dispersion is then written as 
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£n(Pz) = (n + l/2)w c - 2tcosp z s (1) 



where s is the interlayer distance, t the interplane transfer integral, p z the interlayer momentum in z direction, 
n a positive integer, lo c — eH/m* the cyclotron frequency depending on the quasiparticle effective mass m* and 
e the electron's charge (Planck's constant, Boltzmann's constant and the velocity of light are chosen to be unity 
throughout the paper). The Landau quantization of the 2D motion gives rise to the well-known magnetic quantum 
oscillations effects such as de Haas- van Alphen (dHvA) effect - magnetization oscillations - or Shubnikov-de Haas (SdH) 
effect - magnetoresistance oscillations. In quasi-2D organic conductors, the resistivity component in the direction 
' perpendicular to the layers p zz is typically studied [1], when both electrical current and magnetic field are parallel to 
same direction z. 

Recently, J. Wosnitza et al. [2] and M.-S. Nam et al. [3] reported unusually strong oscillations of p zz in the quasi-2D 
organic conductor /3"-(BEDT-TTF)2SF 5 CH 2 CF2S03 in high-magnetic fields at low temperatures. The authors of 
Ref. [2] have called their observation a field-induced metal-insulator transition. M.-S. Nam et al. [3] proposed that 
the insulating-like behavior stems from the presence of a gap in the density of states by pointing out the thermally 
activated behavior of maxima of p zz (or minima of conductivity a zz — p zz ), which occur when the chemical potential 
p, lies between Landau levels. Their interpretation is based on the existence of negligibly small chemical potential 
oscillations with magnetic field owing to the quasi-lD part of Fermi surface which plays the role of an infinite reservoir 
of states. The assumption of a fixed chemical potential is consistent with the information provided by the inverse 
sawtooth wave form of the magnetization oscillations observed in the same compound [2,4], which is a strong indication 
of negligibly small chemical potential oscillations (see e.g. [5] and references therein). 

Some articles [6-10] are devoted to the theoretical study of the oscillations of magnetoconductivity a zz in quasi-2D 
metals. Very recently, P.D. Grigoriev et al. have proposed in a series of papers [7-9] an explanation of other features 
concerning a zz oscillations such as the field-dependent phase shift of the beats and the slow oscillations of a zz , in 
frame of both the semiclassical Boltzmann theory [7] and the quantum transport theory [8,9]. They showed that 
these features observed in layered organic conductors other than /3"-(BEDT-TTF)2SF5CH2CF2S03 originate from 
the pronounced low-dimensionality of energy spectrum (1) in the low-field limit uj c < Airt. It is worth noting that 
the latter condition is not fulfilled in the range of magnetic fields where experiments of Ref. [2,3] were done. The 
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observation of very sharp extrema in the dHvA oscillations [2,4] hints that the 2D limit 4irt <C lo c is most probably 
reached [5]. 

In this paper, we investigate theoretically the magnetoconductivity oscillations with a view to give an explanation to 
the experiments reported by J. Wosnitza et al. [2] and M.-S. Nam et al. [3]. General calculations of SdH oscillations of 
magnetoconductivity o zz in quasi-2D metals are performed in frame of the quantum transport theory for short-range 
impurity scattering. In the second section, the general expression for <j zz is represented in a more convenient form 
valid for ^ » uo c and an arbitrary self-energy. We demonstrate in Section III the existence of a thermally activated 
behavior for the minima of magnetoconductivity in high quality quasi-2D metals in high magnetic fields. This behavior 
is shown to be not related to the presence of a gap in the density of states but rather to the presence of a pseudo-gap 
in the zero-temperature magnetoconductivity function <r zz {p/uj c ). In Section IV, the impurity self-energy for the 
layered quasi-2D metals is investigated in frame of the self-consistent Born approximation. Finally, in Section V, our 
results are discussed in light of the experimental observations in ^"-(BEDT-TTF^SFsCBbCFsSOa [2,3]. The detailed 
description of the studied microscopic model and the corresponding transport calculations are given in Appendix. 



II. THE LONGITUDINAL CONDUCTIVITY a Z z 



In absence of vertex corrections, the longitudinal conductivity a zz in quasi-2D metals under a quantizing magnetic 
field perpendicular to the layers can be expressed as (see Appendix) 
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where summation has to be done over the energy e and the magnetic quantum numbers n of the Landau levels, the 
spin projections a — ±1, and the interlayer momentum p z . Here 3? means the real part, go is the 2D density of states 
in absence of magnetic field, n' F {e) is the derivative of the Fermi-Dirac distribution function 

n'jp(e) = \ , (3) 

FV ' 4Tcosh 2 ( £ /2T)' V ' 

and T the temperature. The limits of the integral over p z are [— ir/s, +tt/s], and those of the integral over e are 
[— oo, +oo]. The functions G R and G A are the retarded and advanced Green's functions defined by 

G<j,n,pz,€ — {^cr,n,p z ,e) ~ _ r T --p i (4) 

where 

£cr,n, Pz = (n + uj c - 2t cosp z s - \i - u^, e H, (5) 

H e is the electron's magnetic moment, and S(e) = iV e is the energy-dependent impurity self-energy. The real part of 
self-energy E is included in chemical potential yu. In our model, S is independent of the quantum numbers n,p Z7 a 
because we consider only point-like impurity scattering. Lastly, v z is the velocity in the z direction 

Vz(Pz) = a f" Vz = 2stsinp z s. (6) 
dp z 

Formula (2) is one half smaller than that used by Grigoriev et al. [9] without derivation. 

In Eq. (2), two terms contribute to a zz : one originating from the product of the retarded and the advanced 
Green's functions (which we call RA), the other from the product of two retarded (or two advanced) Green's functions 
(contribution RR). Further, we choose to calculate the two terms RR and RA separately in order to see clearly how 
they contribute in the final expression: 
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(8) 



In contrast with the case of metals without magnetic field where the contribution RR is absent (see e.g. B.L. Altshuler 
et al. [11]), we show in the third section that the term RR contributes here to the total conductivity, and specifically 
in the 2D limit. 

Now, we rewrite expressions (7) and (8) by performing the summation over the Landau levels n. For convenience, 
we neglect henceforward the spin splitting, which can be easily introduced in the following calculations. To transform 
the sum over Landau levels, we use the Poisson summation formula 

£/(n) = E / duf(u)e-^, (9) 

n=0 l = -oo Ja 

where a is any number between —1 and 0, to obtain for the first contribution RA 
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Changing the variable u to x — £, u ,p z ~ £, and then setting the lower limit of the integral over variable x to — oo since 
we consider the relevant regime \x 3> uj c for SdH oscillations, we get 
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A straigthforward residue calculation for the integral over x yields 
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Collecting positive and negative l 7 and introducing the first order Bessel function J\{x) denned by 



^^sinVe— = *M. 



we rewrite Eq. (12) as 
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where 
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= e 2 . 90 = e 2 ffo (v 2 )ro (14) 
1 

is the longitudinal conductivity in absence of magnetic field, r = 1/2t the imaginary part of the impurity self-energy 
at zero magnetic field, To the mean free time of quasiparticles also at zero magnetic field. As for (Trr, we perform the 
sum over the Landau levels in the same way to get 

ann = ,0 jds (-n' F (e)) ^ £ (-1)' Ji (^) e~ 2 ^ cos (**^) ■ (15) 
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Let's compare both contributions RA and RR. If we consider the limit of zero magnetic field, the contribution <jrr 
vanishes and we recover the usual Drude formula given by 

a?.?. = ora = co- (16) 

We find again that in absence of magnetic field the contribution RA corresponds to the semi-classical contribution 
to transport [11]. The term RR has no semi-classical equivalent. In the 3D limit 4-7rf 3> lo c or in the dirty limit 
27iT e 3> u> c , the contribution RR represents actually a small correction compared to main contribution RA. However, 
at arbitrary magnetic field, the comparison of the two contributions is not trivial, since both terms are expressed as 
alternated series. We demonstrate in next section that the contribution RR becomes quite significant specifically in 
the two-dimensional limit when 47rf <C co c . 



III. MAGNETOCONDUCTIVITY IN THE TWO-DIMENSIONAL LIMIT 



In the 2D limit uu c » 4nt, the two contributions <jra and <jrr (given by Eq. (13) and (15)) get simplified since the 
Bessel function factors can be expanded as J\(x) ~ x/2 for small values of argument x: 
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In this limit, the field oscillating part of a zz becomes independent of the interlayer hopping integral t, which appears 
only in the zero- field normalization factor <jq. This clearly establishes the physical origin of a zz oscillations as the 
quantization of motion in the plane perpendicular to the magnetic field, which affects the interlayer transport through 
energy spectrum (1). 

In order to compare both contributions <jra and <jrr at arbitrary magnetic field such that u c 3> 47rf, summation 
over the integers I is performed for each contribution. Then, the total magnetoconductivity is written as 
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<Jrr{s) = -a 

Uc (cosh(27r^)+cos(2^ I i-)) 

At zero temperature, magnetoconductivity oscillations are then given by 

<Jzz = vra(p) + vrr(p) = o zz {pi). 



(21) 



(22) 



As already noted at the end of Section II, the term <trr contributes negligibly to the total conductivity a zz in the 
dirty limit 27rr e 3> oj c . From Eq. (21) we can see that at zero temperature <Jrr gives always a negative contribution 
to <r zz {p) when /j, ~ nu c , and a positive contribution when jj, <~ (n + l/2)w c . 

When uo c 3> 2-7rr e , uo c 3> 47rt and \i ^ (n+ l/2)uo c , both contributions RA and RR are of the same order so that the 
total magnetoconductivity at zero temperature vanishes periodically, a zz (p.) — > 0. Such a limit indicates the existence 
of a particular behavior for the form of SdH oscillations at low temperatures, which is the subject of the following 
development. 
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Let's consider that the oscillations of self-energy are negligibly small, i.e. T £ w Tq. This assumption is not a 
necessary ingredient for our theory. We only use it to study the behavior of a zz in high magnetic fields on a more 
quantitative level. However, the relation T £ w Tq is proved to be correct in particular in the self-consistent Born 
approximation in presence of an internal reservoir of states in the conditions of experiment [3] (see next Section) . 

The numerical analysis of the function a zz (e) reveals two regimes depending on the ratio a = 2nT n /uj c for the 
shape of oscillations: 

• A low-field regime u c < 2nT (or a > 1), where the contribution RR represents a correction of the RA 
contribution for the whole range of energy e, 

• A high-field regime defined by uo c > 2-7rro (or a < 1), where a pseudo-gap centered on the energy e = nuj c exists 
in the function a zz {e). 

In Fig. 1 we represent the function a zz (e) for different values of magnetic fields, using the experimental values given 
in Ref. [3], namely T =0.22 meV and m* — 1.96m e (here m e is the electron mass). As shown in Fig. 1, the opening 
of a pseudo-gap occurs for a magnetic field of about 20 T, corresponding to a parameter a ~ 1. The width of the 
pseudo-gap clearly increases with the magnetic field. 

Such a high-field behavior can also be found from analytical considerations. Indeed, the presence of a pseudo-gap 
in CTzz(e) is clear if we write the total conductivity as 



For e <~ noj c and a < 1, magnetoconductivity a zz (e) becomes very small over a range in energy, which gets more and 
more wide with a decreasing a. 

Such small values for a zz when a < 1 lead to huge maxima of magnetoresistance p zz at low temperatures when 
n/uo c ~ n. To illustrate the behavior at finite temperature, we have calculated numerically a zz using Eq. (19)-(21) 
still taking experimental values [3]. We estimated the chemical potential at /i~ 11.8 mcv (for the numerical study, 
we have explicitly assumed that the chemical potential oscillations are negligibly small; effects of chemical potential 
oscillations are evoked in the dicussion). In Fig. 2, the general form of p zz = <j~} oscillations for negligibly small 
self-energy oscillations is represented in the range of magnetic fields 20 - 60 T. In this plot, we have chosen exactly the 
same temperatures as in Fig. 1 of Rcf. [3] in order to compare more easily the shapes of magnetoresistance oscillations. 

As in Ref. [3], we have also studied numerically the temperature dependence of the resistivity maxima. The 
presence of a pseudo-gap for high-magnetic fields gives rise at finite temperature to a thermally activated behavior 
of conductivity minima at fj, ~ nuj c . Indeed, according to Eq. (19) and when fi ~ nw c (then the minima of a zz (e) 
are in correspondence with the maximum of —n' F (e)), the main means of conduction is the thermal excitation of 
quasiparticles at the edges of the pseudo-gap. Figure 3 shows \na zz versus 1/T for the integer values n/w c = n =4, 
5, 6, 7 and 8. In this plot, a linear dependence is clearly found for the interval of temperatures 2 < T < 4 K (or 
0.25 < T~ x < 0.5 K _1 ), pointing out the activated behavior of a zz minima. At very low temperatures T < 2 K 
(or T^ 1 > 0.5 K _1 ), a zz saturates (Fig. 3): quasiparticles whose energy lies in the pseudo-gap mainly contribute to 
conduction. 



For the numerical study of preceding section, we have assumed that the oscillations of the self-energy arc negligibly 
small in the 2D limit uo c 3> 47rf. The purpose of this section is to justify such a choice in the organic compound 
/3"-(BEDT-TTF)2SF5CH2CF 2 S03. The impurity self-energy is derived in the so-called self-consistent Born approxi- 
mation. In quasi-2D metals, we usually have to take into account of the presence of an internal reservoir of states [6,4]. 
For short-range impurity scatterers, the imaginary part T of the self-energy is directly proportional to the total density 
of states g tot , which is taken as the sum of the quasi-2D density of states g and of the field-independent reservoir's 
density of states g R . We assume that the same relation of proportionality holds both in presence and in absence of 
magnetic field and thus 
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IV. THE IMPURITY SELF-ENERGY 




(24) 
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where R = g /go is a parameter measuring the strength of the reservoir, g(e) the density of states of the quasi-2D 
metal in presence of a quantizing magnetic field, and go = go/ s the quasi-2D density of states at zero magnetic field. 
For simplicity, we neglect the dependence of the reservoir's density of states g R on energy. The quasi-2D density of 
states g{e) depends on the impurity self-energy through the expression 



(25) 



which for /i>w c can also be written (without spin-splitting) after application of the Poisson summation formula 
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Then, after summation over momentum p z we obtain [12] 
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where Jo is the zero order Bessel function. Combination of Eq. (24) and Eq. (27) yields a complicated self-consistent 
equation to solve for the impurity self-energy T e . 

For 47rf > cj c , it is reasonable to keep only the first term I = 1 in the sum of Eq. (27). As a result, the impurity self- 
energy oscillates weakly with H. In absence of reservoir (R — 0), this field dependence has to be taken into account 
in the expressions (13) and (15) of the longitudinal conductivity a zz . This regime has already been investigated by 
Grigoriev et al. [7-9] who showed that the field-dependent phase shift of the beats and the slow oscillations of a zz are 
consequently general features of the SdH oscillations in low-dimensional metals. 

In general, we have to take into account many terms in the sum of Eq. (27). By summing over the integers / before 
summing over momentum p z in Eq. (26), we can write the quasi-2D density of states g(e) in a more convenient form 
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so that the self-consistent equation to solve in quasi-2D metals becomes 
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For uj c <C 27rFo, the Landau levels are strongly broadened by impurity scattering and the oscillations of impurity 
self-energy scattering are negligibly small; this is not the interesting regime. Expression (29) enables to study the 
impurity self-energy at arbitrary magnetic fields. The self-consistent equation to solve in the 2D limit uu c 3> 47rf takes 
the form 
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sinh ( 27T— 



To 1 + R 1 + R 



(30) 



This equation in absence of reservoir (R = 0) is similar to that obtained many years ago by T. Ando for a 2D 
metal [13]. For magnetic fields lu c > ttTo, it leads formally to a zero value for the self-energy T e when energy s lies 
between two adjacent Landau levels [13]. But for T £ = 0, one can not fulfill the summation in (26) (taken with t = 
for the 2D metal) . This contradiction demonstrates that the self-consistent Born approximation is inappopriate in the 
high magnetic field region lu c > ttTo- To avoid this unphysical result obtained for a short-range impurity scattering, 
M.E. Raikh and T.V. Shahbazyan [14] (see also Ref. [15]) have discussed the 2D electron system in a random potential 
with a finite correlation radius R c . It was shown that the self-consistent Born approximation is justified only in the low 
field region when the magnetic length ^ < i! c . At high magnetic fields such that R c > Ifj, the higher order diagrams 
contributions become important and have also to be taken into account in the impurity self-energy calculations. As 
a result, the density of states is proved to be finite at any energy and for arbitrary magnetic fields. 
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However, for R ^ 0, we note that the unphysical cut-off for the impurity self-energy does not exist. This can be 
shown by rewriting Eq. (30) as 



This equation admits solutions for the whole range of energy e and for any magnetic fields when T e > TqR/(1 + R). 
Thus, it is always possible to find a finite impurity self-energy at arbitrary energy and magnetic fields. For this reason, 
we consider in the present paper that in the 2D limit the oscillations of T e are given by Eq. (31). The validity of the 
self-consistent Born approximation in the 2D limit and in presence of a reservoir is addressed to future work. 

We plot in Fig. 4 the oscillations of r e for different magnetic fields H =20, 40 and 60 T, and for a reservoir 
parameter R = 5 (the choice of such a value was suggested by Ref. [4]), using the same parameters as in Fig. 1. As 
expected, the oscillations amplitude increases with the magnetic field (see Fig. 4). However, for R = 5, oscillations 
of self-energy r e are relatively small compared to To even at magnetic fields of about 60 T. In absence of other 
informations about the importance of the reservoir's influence on self-energy oscillations, we consequently assume at 
this level that oscillations of T e are small in experiments of Ref. [2,3]. 



Numerical treatment of Eq. (19)-(21) allows comparison with experiments [2,3]. Theoretical calculations with 
negligibly small self-energy oscillations reproduce well the appearance of huge oscillations in the same conditions of 
temperature and magnetic fields as in experiments [2,3]. As demonstrated in third section, the thermally activated 
behavior of p zz maxima is theoretically established in correspondence with the observations of Ref. [3]. Furthermore, 
the minimum value of magnetic field to have thermal activation, i.e. 27rFo ~ u) c , is in quantitative agreement with 
experiment (in the conditions of experiment [3], it corresponds to a magnetic field of about 20 T). 

Nevertheless, we notice that there exist some discrepancies between our calculations and experiment [2,3] especially 
at very low temperatures. Indeed, our numerical calculations do not qualitatively reproduce the behavior at very low 
temperatures of the resistivity minima, which are increasing with the magnetic field (compare Fig. 2 of the present 
paper and Fig. 1 of Ref. [3]). Quantitative departures from experiment are also found for the resistivity maxima at 
very low temperatures. A reason for the latter may be the consideration of negligibly small self-energy oscillations in 
our calculations. Indeed, although the presence of self-energy oscillations does not suppress the thermally activated 
behavior of p zz maxima, it induces quantitative modifications for the pseudo-gap value and thus for the resistivity 
maxima. However, this can not qualitatively explain the field dependence of the resistivity minima at very low 
temperatures. 

A possible limitation of our model to describe the experimental signal [2,3] is the consideration of point-like impurity 
scattering, which enables to neglect the vertex corrections. Indeed, starting from a microscopic model with the energy 
spectrum (1) (i.e. a 3D model), we find that in presence of point-like impurity scattering the behavior of the magnetic 
quantum oscillations of a zz is fully dictated for magnetic fields oj c 3> 47rf by the 2D in-plane motion. Under the same 
condition and in a parallel direction, the equation for the density of states of the quasi-2D metal investigated with the 
self-consistent Born approximation takes also a 2D form. On the other hand, in the pure 2D case, the treatment of the 
impurity effects on the magnetic quantum oscillations by a short-range impurity scattering is known to be no more 
valid at high magnetic fields to describe accurately the density of states. According to Ref. [14], when typically the 
magnetic length Ih < R Cl it becomes necessary to include also the effects of the impurity scattering on finite range R c 
in the calculations of the density of states. Probably, specific contributions due to the finite-range impurity scattering 
have also to be taken into account in the calculations of the magnetoconductivity oscillations in high magnetic fields. 
In our opinion, such an open question might explain why our present calculations do not seize the complete correct 
behavior of magnetic quantum oscillations of p zz in the high field regime. 

So far, we did not discuss the possible effects of the chemical potential oscillations on the magnetoconductivity 
oscillations. As declared in introduction, the dHvA signal [4] gives strong support for negligibly small chemical 
potential oscillations in the organic compound /3"-(BEDT-TTF)2SF5CH2CF2S03, what we use at the final stage 
to perform the numerical calculations. The mechanism responsible for fixing the chemical potential is not clearly 
established. Actually, this question is directly related to the theoretical problem just evoked above of the density of 
states of a 2D metal in high quantizing magnetic fields. Nevertheless, it is interesting to note that the value of the 
pseudo-gap appearing in the zero-temperature magnetic field-dependence of a zz depends strongly on the presence of 
chemical potential oscillations. In third section, it was shown that in the limit uj c — > oo the magnetoconductivity 
&zz{p) vanishes everywhere but at p, = (n + l/2)u) c . In a pure metal in absence of reservoir, the chemical potential is 




(31) 



V. DISCUSSION 
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pinned to the Landau levels, i.e. fj, = (n + l/2)tu c when varying the magnetic field. Therefore a zz does not vanish, 
except when the chemical potential jumps. Thus, in the pure (unrealistic) case, the pseudo-gap in a zz is reduced to 
a point and there is no thermally activated behavior for the conductivity minima. In a real system, the chemical 
potential oscillations at zero temperature lead to a smaller pseudo-gap value than that obtained at a constant chemical 
potential. Although no full quantitative agreement between our numerical calculations and experiment has been found 
for the pseudo-gap value, the observation [3] of a well developed thermally activated regime with a relatively important 
pseudo-gap value seems accordingly consistent with the existence of negligibly small chemical potential oscillations. 

We point out that the observation of a drastic change in the temperature and magnetic field behavior of <j zz 
oscillations is a unique feature of the 2D limit with energy spectrum (1). Since the switch between the two regimes 
occurs for w c ~ 27rr , this necessarily implies that T > 2t. Such a condition is exactly expected to be reached in the 
organic compound /3"-(BEDT-TTF)2SF 5 CH 2 CF2S03 [2]. It is worth noting that the same condition is proposed for 
having incoherent intcrlayer transport [2]. In that case, the existence of an interlayer dispersion in the energy spectrum 
(1), or in other words the Fermi liquid description, is questioned. It was with such an idea that authors of Ref. [2] 
attempted to extract a background magnetoresistance superposed with SdH oscillations. In a recent experiment [16] 
performed with the same organic compound as concerned here, another arguments in favor of incoherent transport 
have been proposed to interpret observations of the angular dependence of p zz in high magnetic fields. However, we 
notice that the coherent magnetotransport in quasi-2D compounds is usually discussed in frame of the Boltzmann 
transport theory. Here, we clearly find important discrepancies between the results obtained with the quantum 
transport theory and those obtained with the Boltzmann transport theory especially in high magnetic fields io c ^> Ant 
and u c > rV Therefore, theoretical studies developed exclusively in frame of the quantum transport theory are needed 
to better understand the possible unconventional quantum features of coherent magnetotransport in low-dimensional 
systems. As for the possible absence of the quasi-particle description at t < t^ 1 in the quasi-2D metals [16], one 
can notice that there is no doubt concerning the quasi-particle concepts in the pure 2D case where t = 0, and the 
only criterion for this rests on the Fermi energy: ep > Tq 1 . It seems that the possibility of an interlayer motion at 
t docs not spoil the quasi-particle description because the in-plane scattering just leads to the undetermination 
5k <~ (vfto) of the in-plane momentum (vf is the velocity at the Fermi surface) but not to the indetermination of 
the interplane quasi-momentum k z . This explains in our opinion the observations of the coherence in the quasi-2D 
organic conductor k-(BEDT-TTF) 2 Cu(NCS) 2 reported at t^ 1 - 6t in Ref. [17]. 

Finally, the importance of the contribution (jrr in the 2D limit and in high-quality layered conductors demonstrated 
in Section III prevents from separating the magnetoconductivity o zz into a field oscillating part and a steady part in 
the manner used by the authors of Ref. [2]. The relation of proportionality between the derivative of dHvA signal 
and SdH signal assumed in Ref. [2] breaks down in high magnetic fields and in low-dimensional systems. This already 
occurs for w c < Ant due to small self-energy oscillations in absence of reservoir [7-9]. Furthermore, even in presence 
of negligibly small self-energy oscillations, the relation of proportionality may not be valid. In this latter case the 
integration over energy e in (17) can be performed analytically and yields for the oscillating part of the contribution 
RA in the 2D limit: 

<Jra - <r = 2a Q J2(-l) 1 e- 2nl ^ cos (^irl-^ Jde(-n' F (e)) cos (^nl^j (32) 

+ oo , 

= 2^ (-l)V 2 ^cos(2^ 
i=i V Wc 



sinh A; 



(33) 



where we have introduced 



2ir 2 lT , , 

A, = . (34) 



In quasi-2D metals, the magnetization oscillations M are given in the 2D limit by [12] 



M = go ^--^ e -sin [2nl- . (35) 



H tt I sinh Xi \ w, 

From Eq. (33) and (35), we can derive (in the case of negligibly small chemical potential oscillations) a relation valid 
at any temperature between the derivative of M and the oscillating part of the contribution RA by keeping only the 
rapidly oscillating terms: 

g ra - cr H 2 dM 

~ 2 tef- y Sb > 

cr g fi 2 dH 
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In the dirty limit u c -C 27rro, it has been shown in section II that ctrr contributes negligibly to total conductivity, 
i.e. a zz i=a <7ra, so that the SdH oscillations are approximatively proportional to H 2 dM/dH , exactly as in 3D metals. 
When u c > 2nT , the contribution a^R becomes quite important and there is then no more simple correspondence 
between dHvA and SdH signals. 



VI. CONCLUSION 



We have investigated the oscillations of the longitudinal magnetoconductivity in quasi-2D organic metals for im- 
purity scattering on short-range in a perpendicular quantizing magnetic field. Huge magnctoresistance oscillations 
with unconventional magnetic field and temperature dependences have been found in high magnetic fields and for low 
temperatures in layered conductors whose interlayer hopping integral Airt <C lo c . In particular, we showed that the acti- 
vated behavior of resistivity maxima recently observed in the organic conductor /3"-(BEDT-TTF)2SF 5 CH2CF 2 S03 [3] 
is a general feature of Shubnikov-de Haas oscillations in the 2D limit oj c » Ant in good samples and/or in very high 
magnetic fields lo c > 27iTo. The origin of such a behavior is not related to the presence of a gap in the density of 
states but rather to the presence of a pseudo-gap in the zero temperature magnetoconductivity function a zz {ji/u> c ). 



APPENDIX: CALCULATION OF o Z z 



The longitudinal conductivity a zz characterizes the appearence of a current density in the direction z in response 
to an external electrical field E z 

jz=cr zz E z . (Al) 

With the objective to find the longitudinal conductivity, an alternative monochromatic vector potential A z depending 
on the Matsubara variable r is formally introduced in the form 

A z {r,j,r) = A z {k,k z ,ujy % - F+ik ^-^ r (A2) 

and such that 

E, - -^i. (A3) 

Here, j is an integer which indexes the layer and r — (x,y) is the position in the layer. The relation between E z and 
A z takes a simpler form in Fourier space 

E z (k,k z ,uj) = iu>A z (k,k z ,u)). (A4) 
Then, the relation between the formal field A z and the current density in Fourier representation is 

j z (k,k z ,w) = J z (r,J,r)e^- p ~ tk ^ + ^d 2 rdr = a zz icoA z (k, k z ,co) (A5) 



1. Hamiltonian 



We consider the same model Hamiltonian as that used in the paper of R.A. Klcmm et al. [18], in which the electrons 
are free to propagate within each layer, and tunnel between the layers (in the vertical direction z). The Hamiltonian 
is constructed in the second quantization language, where the field operator ^ j,o-(r) corresponding to an electron with 
spin a at position r in the jth layer is defined by 

^ j Ar) = ^=Y j e it %Ak), (A6) 
% 

where aj i<T (fc) is the annihilation operator and 5* is the area of a layer (S = 1 after). The Hamiltonian we shall consider 
is of the form 
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H — Ho + Ht + H.% 



(A7) 



where Ho is the two-dimensional Hamiltonian for an electron in a magnetic field H = Hz associated with a vector 
potential A 



H = sJ d 2 r]T [*t CT (r) (-»Vr- eA) - ^ e a ■ H^j *^(r) 



Ht is the gauge-invariant tunneling Hamiltonian, which models the coupling between the layers : 
H T = -ts J d 2 r^*] ff (f)^ +l!(T (r)exp (^ le ' A z (f,j,T)dz^j +H.C. 



(A8) 



(A9) 



Here the expression H.C. means the Hermitian conjugate. Finally, Hi mp describes the coupling of electrons with the 
impurities in the layers (we neglect the scattering that might occur in the tunneling process) 



H imp = s J (frj^^l^u^j^r), 



(A10) 



where Uj (r) is the short-range scattering potential. 

In our model, the longitudinal current density j z originates only from the tunneling process, and is then given by 



3z{r,j,r) 



SHi 



5A z {f,j,T) 

The derivation of the Hamiltonian Ht yields for a small external field A z 

Jz(r,j,r) = -iets 2 J2 [^>(0^+i,a(r)exp(-ise^(f,j,r)) - ff.C, 

<7 
<7 



(All) 



(A12) 



(A13) 



2. Green's functions 

To proceed further, we shall use the usual temperature Green's functions defined by [19] 

Gl^{r uTl ,r 2 ,r 2 ) = - (T^^irun)^^,^)) (A14) 

where the averages are evaluated in the grand canonical ensemble, and (r\, t{) are the Heisenbcrg field operators 
depending on the Matsubara variable t\. Therefore, writing the averaged current density in the jth layer in terms of 
the Green's function of the system, we obtain 

h t) = (j z (r, j, t)) = jf (r, t) + jf (r, r) (A15) 

where 

jf (r,r) = -iets 2 J2 (tfZijfr t,P,t) GfJ +1 (f,T,r,r)) (A16) 

is the paramagnetic part of the current density and 

jf (r, r) = -e 2 ts 3 A z (r, j, r) £ (GfJ +1 (r, r, f, r) + Gf*,. (r, r, f, r)) (A17) 

cr 

is the diamagnetic part of the current density. 
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The equation of motion for the Green's function in the presence of impurities, magnetic field, and the vector 
potential A z (f,j, r) is, for a given spin a 

- Uj (r)} G a ^(r,T^,T>) = s^SjtSir- 7)5(t - r') (A18) 

where the operators Aj± are defined by 

A J± vEv(r» = f J± i, CT (r,r). (A19) 

The average of Eq. (A18) over impurities configuration is done using ('Uj(r)^^^(f , ')) = UQS~ 1 5jiS(r — f ) (here w is the 
amplitude of impurity potential). Keeping for the averaged Green's function (G) the same notation G as before, we 
get for the equation of motion for the averaged Green's function 

4? ~ i (-^ ~ eA ) 2 + 1 + &i-e iesA °) + M + WH+ 

+Z}G (7J j(r,Ty,T , )= S - 1 S 3l 5(r-?)5(T-T') (A20) 



where E is the impurity self-energy. 

i's functic 

Gv,j,l(r, r, f, t<) = G° ;j _ ; (r, f, r r') + G« ; (r, r, r<, r') (A21) 
in 

It is more convenient to write G° in the Landau representation : 



Next, we write the Green's function Gj.i in the form 



where G?, is the Green's function in the absence of A z and Gy) is the first correction which is linear in A z . 



X 



t/s 

^ * n (a: - p^) * n (a/ - p^) G°,„ ;pz (w„) (A22) 



where 



and 



^ w = i,-u + %) (A23) 



= -— T7T-— ^= exp l~]ff„[^l. (A24) 



Here i? ra is the Hermitte's polynom of degree n, Z# = 1/y/eH is the magnetic length, 

fa,n, P , = + w c - 2t cos p z s - \i - a/j, e H, (A25) 

w„ = 7tT(2v + 1) are Matsubara frequencies and r(w„) is the imaginary part of the impurity self-energy E. The real 
part of E is included in chemical potential. Here we consider that E is independent of the quantum numbers, which 
is valid for short-range impurity scattering. 
As for G^\ it is given by 

G%{r,r,?,r') = sY, ff^r"dr"Gl tj _ k (f, r", r - r")V k {f , ^" ~ ^ ( A26 ) 

where 14 is the potential 

V fe (f,r) = ietsA(r,fc,r) (A fe+ - A fe _) (A27) 

which is taken as a perturbation. 
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3. Diamagnetic current 

At first order in the vector potential A z , the diamagnetic current is thus 

jf (r,r) - -e 2 ts 3 A z (r,3,T)Y, (G^r, f, 0) + G° 4 (r,f,0)) (A28) 

From Eq. (A22), we have straightforwardly 

Gl^iW) = f c c t]T f " |£e*"W-'> £G°,„, p >„) (A29) 



where 



9 V = ^ = /tr [*-(*- 2 (A30) 



2 2tt^ ./ 2tt 



is the degeneracy of Landau level per spin state. Then, the diamagnetic current is 

jf (r,r) = -e 2 t S 3 A(f,j,r) 3oWc E | ^ cos(p z a)r^ G°^ ( W „). (A31) 

Next, we transform the summation over Matsubara frequencies into an integration, and perform the analytical 
continuation [20] 

T^GwW = JZ / def e (G* nj>MiB - G^ p J (A32) 

^ n J -co 



where 



f e = tanh ±. = \-2n F (e). (A33) 



Then, rewriting the expression for the diamagnetic current density by integrating by parts over the variable p z , and 
making the Fourier transformation, we find 

j D (k,k z ,u J )=A z (k,k z ,u J yA c Y: J f^2(p,) / ^J,{(Gl ntP J 2 -(G^ p J 2 }. (A34) 



4. Paramagnetic current 

The paramagnetic part of the averaged current density is given by 

jf (f, r) = - le ts 2 ]T [g« +1j (r, r, r, r) - G« J+1 (r, r, f, r)] . (A35) 

We are looking for the Fourier transform of the current density 

jf(k, k z ,co) = s JJSrdre^-^Y^e-^jfir, r). (A36) 

The summation over j yields 

*E e " iMj ( G ?+U - G S + i) = / - ^ {P ^ )S )G { Z^. (A37) 

3 

By definition, we have 
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G £U( r >' r >) =^ e -^^-G«(r-',r,r,r). (A38) 
Using (A26) and (A27), we obtain 

G £x( r >> r >) = - ffd 2 r"dT"Gl(ry\T-TnevM)Mr 7 \p z -p' z ,TnGl,(7',r,T" - r). (A39) 

Making the inverse Fourier series transformation with respect r and the Fourier series transformation with respect to 
r, we have 

d 2 rdTe^ T '^G^\ p , (r, r, r, r) = -|u; c T £ G°, p , (u + )ev z (p' z )A z (k,p z - (A40) 

where w± = w„ ± w/2. Collecting (A36), (A37) and (A40), we finally obtain 

j p (fc, k z , w) = %A Z (k, k z , u)e 2 ts 2 f^cE/^ S - ^ > s ) v z (p z - k z ) x 

n.cr 

xr^G° Wj ( U+ )G Wi . l| ( W .). (A41) 

We are interested in the limit £^ <C because E z is assumed to be uniform. Then, 

j p (k,k z ,u) = -A(fc,fc„c)e 2 S |c c E / |^( Pz )T^G° n , Pi K)G^ !P! (^). (A42) 

Next, after transforming the summation over Matsubara frequencies into an integration over energy e and taking the 
analytic continuation [20], we get 



efe 
47ri 



G pz,e ife fe-u))G PzS _ UJ -\-G p ^ e f £ - u G p ^ e _ UJ G pzE f e G pze _ ul \. (A43) 



5. Static conductivity 



At zero order of u>, the paramagnetic current density arising from the term RR is cancelled with the diamagnetic 
current density (see Eq. (A34)) 



& M +i& = 0. (A44) 

At first order of u, a nonzero total current density arises from the terms RA + RR of the paramagnetic current. Using 
Eq. (A5), we find the static conductivity 

a„ =e 2 af c c ]T f f^ 2 fe) / g {f'^G^ + 2fJSt (<£,.<#,.)) > ( A45 ) 

Tl,<T 

where the prime means the derivative with respect to the energy e. Integrating by parts over the energy in the term 
i?_R and using (A33), we obtain Eq. (2). 
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FIG. 1. Presence of a pseudo-gap in the function a zz (e) for high magnetic fields. The opening of a pseudo-gap occurs for 20 
T, i.e. for a parameter a — 2-kF /lu c « 1. Thus, the other magnetic fields H — 10, 30 and 40 T correspond respectively to a = 
2, 2/3 and 0.5. 
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FIG. 2. Oscillations of the magnetoresistance p zz /po = <To/cr zz as a function of magnetic field for different temperatures 
(from the top, 0.59, 0.94, 1.48, 1.58, 1.91, 2.18, 2.68, 3.03, 3.38, 3.80, and 4.00 K). 
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FIG. 3. Plot of \na zz /ao versus 1/T (T in Kelvins) for the integers n = h/lo c (maxima of resistivity). 
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